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Abstract 

Background: The relationship between serum biomarkers and clinical expressions of COPD is limited. We planned 
to further describe this association using markers of inflammation and injury and repair. 

Methods: We studied lung function, comorbidities, exercise tolerance, BODE index, and quality of life in 253 COPD 
patients and recorded mortality over three years. Serum levels of Interleukins 6,8 and16, tumor necrosis factor alpha 
(TNF a) [inflammatory panel], vascular endothelial growth factor (VEGF), and matrix metalloproteinase 9 (MMP-9) 
[injury and repair panel] and pulmonary and activation-regulated chemokine (PARC/CCL-18) and monocyte 
chemotactic protein 1 (MCP-1/CCL2) [chemoattractant panel] were measured. We related the pattern of the 
biomarker levels to minimal clinically important differences (MCID) using a novel visualization method [Observed 
Clinical Association Results (OSCAR) plot]. 

Results: Levels of the inflammatory markers IL-6, TNF a were higher and those of injury and repair lower (p < 0.01) 
with more advanced disease (GOLD 1 vs. 4). Using the OSCAR plot, we found that patients in the highest quartile 
of inflammatory and lowest quartile of injury and repair biomarkers level were more clinically compromised and 
had higher mortality (p<0.05). 

Conclusions: In COPD, serum biomarkers of inflammation and repair are distinctly associated with important 
clinical parameters and survival. 
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Introduction 

Chronic obstructive pulmonary disease (COPD) is a 
complex inflammatory lung disease with systemic reper- 
cussions that impact on patient's performance and sur- 
vival [1]. The clinical presentation, disease severity and 
progression are quite heterogeneous and maybe the 
result of a diverse pathogenic processes that involves 
abnormalities in different pathogenic pathways (protease- 
antiprotease balance, oxidative distress, altered immune 
system and/or apoptotic control) These pathways and 
other mechanisms may act distinctly or in concert in 
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individual patients and be responsible for the different 
phenotypic expressions of the disease [2]. 

COPD is thought to be intimately linked to inflamma- 
tion, documented locally and systemically [3] and as a 
consequence, attention has been centered on the level of 
inflammatory markers and their relation to clinical and 
physiological measurements [4-6]. 

Among them, serum biomarkers have been increas- 
ingly described in cross sectional or short interventional 
studies and associated to meaningful clinical outcomes 
[7,8]. Plasma C- reactive protein (CRP) level has been 
found to be associated to disease severity, quality of life, 
exercise capacity, response to treatment and mortality 
[9,10]. Fibrinogen has also been well studied and asso- 
ciated to survival, risk of exacerbation and poor clinical 
outcome [11,12]. Other studied biomarkers include Des- 
mosine and Isodesmosine ( markers of degradation of 
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mature elastic fibers) [13] serum amyloid A [14], pro- 
adrenomedullin [15], procalcitonin and CXCL-10 (during 
exacerbation) [16,17], surfactant protein D [18], serum 
PARC/CCL-18 [8], CC-16 [19]and fibronectin [5]. Not- 
withstanding, there is limited association to disease sever- 
ity (except for CRP) and to other clinical outcomes [20]. 

Using high throughput proteomics, we showed that 
the serum level of 24 out of 147 analytes separated 
patients with COPD from smokers and non-smoker 
controls [7]. We also associated the level of selected 
biomarkers with lung function, functional capacity, 
exacerbation rate and the BODE index. 

In the present manuscript, we reasoned that the 
phenotypic expressions of COPD may be the product 
of the balance of different biological pathways such as 
inflammation and injury and repair and that a com- 
posite panel of biomarkers expressing different mech- 
anistic pathways could be useful in expressing the 
complex nature of COPD. We also planned to further 
study the association of serum biomarkers with out- 
comes applying more stringent criteria that involve the 
use of minimal clinically important differences in out- 
comes and survival. 

To test these hypotheses we collected baseline serum 
samples from a large cohort of patients with COPD in 
2 different centers. We selected 8 analytes from the 
original 24 biomarker panel that expressed different bio- 
logical pathways and correlated them to clinical out- 
comes using a novel color coded method that relates the 
level of the biomarker to clinical expressions (beneficial 
or not) of the disease and possibly represent different 
pathobiological functions. 

Materials and methods 

This is a prospective cohort study of 253 COPD patients 
representing all stages of disease severity as defined by 
GOLD [21] recruited from two BODE [1] cohort sites 
(St Elizabeth's Medical Center, Boston USA and Hospital 
Universitario Nuestra Senora de Candelaria, Tenerife 
Spain). The Institutional Review Board approved the 
study at both institutions and participants signed the 
informed consent. 

Patients were included if they fulfilled criteria for 
COPD [22] ( post bronchodilator FEVj/FVC < 0.7 and > 
10 pack/year smoking history) and were in stable con- 
dition (no exacerbation for at least 3 months). Exclu- 
sion criteria included myocardial infarction, within 
the past six months, angina, congestive heart failure, 
malignancy, hepatic cirrhosis, end-stage renal disease, 
rheumatoid arthritis, orthopedic condition precluding 
performance of a walking test neurological or psych- 
iatry illness that interfere with participation in the 
study, tuberculosis or other infection that could affect 
biomarkers measurement. 



Participants completed a medical history, physical 
examination, questionnaires regarding anthropometrics, 
respiratory symptoms, co-morbid conditions (Charlson 
Index [23], current medications including inhalers, sys- 
temic corticosteroids, xanthenes, statins, antihyperten- 
sive, exacerbation requiring intervention [24], dyspnea 
(mMRC) [25], and health status (St George's Respiratory 
Questionnaire, SGRQ) [26]. Lung function and two six 
minute walk distance (6MWD) tests were measured 
following the ATS standards and guidelines [22,27]. 
COPD was categorized by GOLD stages [21] and the 
BODE index [1]. Survival was measured using medical 
charts, contacting the patient, family members, primary 
care physician or reviewing a mortality death index 
yearly. 

Biomarker panel 

Based on our previous work [7] and potential disease path- 
ways [28], we selected 8 biomarkers that also showed 
strong association with COPD. They were grouped as 
Inflammatory (IB): interleukins 6, 8, andl6 (IL-6, IL-8, 
IL-16) and tumor necrosis factor alpha (TNF a); Injury and 
Repair Biomarkers (IRB): vascular endothelial growth factor 
(VEGF) and matrix metalloproteinase nine (MMP-9) and 
primary chemoattractants biomarkers (CB): monocyte 
chemoattractant protein-1 (MCP-1/CCL2) and pulmonary 
activation regulated protein (PARC/CCL-18). Blood sam- 
ples were drawn, centrifuged and the serum frozen at -80 
degrees centigrade. The levels were determined using indi- 
vidual ELISA tests performed by a blinded group of 
researchers at the Universidad de Navarra in Pamplona, 
Spain The lower limit of detection for each marker is as fol- 
lows: (IL-6: 0.70 pg/ml, IL-8: 3.50 pg/ml, IL-16: 6.20 pg/ml, 
TNF a,: 1.6 pg/ml, MCP-1: 5.0 pg/ml, PARC: 10 pg/ml, 
VEGF: 5 pg/ml, MMP-9: 0.156 ng/ml). 

Data analysis 

Continuous variables were expressed as means (SD) and 
medians [IQR]; categorical variables were reported as 
proportions. For group comparisons, chi-square was 
used for categorical and analysis of variance (ANOVA) 
or covariance (Kruskal-Wallis test) for continuous vari- 
ables. Statistical significant differences were observed at 
the level of p < 0.05. Associations among biomarkers 
were explored using Spearman test. 

The relationship between biomarker levels and clinical 
outcomes was explored using a novel heat map expres- 
sion. The patients were divided into quartiles by their 
serum biomarker level. Subjects belonging to the upper 
(75* percentile) and lower (25^^ percentile) quartiles by 
their specific biomarker level were compared in terms of 
clinical characteristics and outcomes, including the 
FEVi, inspiratory to total lung capacity ratio (IC/TLC), 
DLCO, dyspnea, health status, 6MWD and BODE index. 
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We assessed the importance of differences in clinical 
outcomes using known minimal clinically important 
differences (MCID) for the variables where it is estab- 
lished [29]: FEVi(100 ml), MMRC (1 point) [30], SGRQ 
(- 4 points), BODE( 1 point) [22], 6MWD (37-71 m) 
and IC/TLC (<0.25) [31]. Statistical comparisons were 
used to test differences for the other variables with un- 
known MCID, including death over 3 years. 

The heatmap colors express the relationship between 
upper or lower quartile biomarker level and the clinical 
outcome (ObServed Clinical Association Results or 
OSCAR plot) using the following color-coding: Differ- 
ences associated with good clinical outcomes are shown 
in green, poor outcomes in red and no differences in 
yellow. 

To determine the relative predictive contribution of 
the biomarkers with mortality as the outcome we com- 
pleted C statistics using only the biomarker panel or as- 
sociating it to the BODE a well-validated predictor of 
mortality in COPD. 

Results 

The baseline characteristics of the cohort are shown in 
Table 1. The patients represented all stages of COPD 
severity (GOLD stage I: 7%, II: 29%, III: 42% and IV: 
23%), 64% were men, mean age 65 (9) years, with im- 
portant hyperinflation (residual lung volume 212 (64) % 
predicted value). Despite severe lung disease, the group 
had a preserved BMI (mean: 27 kg/m^), nutrition status 
(FFMI mean: 19.62), exercise capacity (6MWD mean: 
422 m), reduced comorbidities (Charlson Index mean: 
0.45) and relatively low BODE index score (mean: 3 
points). 

Biomarker values changed significantly with rising dis- 
ease severity. In general, inflammatory and chemo- 
attractant biomarkers increased with worsened disease 
severity despite more frequent use of systemic corticos- 
teroids while biomarkers of injury and repair decreased 
(Table 1 and Figure lA-C). Other medication use (inha- 
lers, statins, antihypertensive medications) was similar 
between patient's quartiles. 

Biomarkers associations 

The relationship between individual biomarkers is 
shown in Table 2. The associations were significant 
among some of the markers. Inflammatory-pathway bio- 
markers IL-6, IL-8, IL-16 and TNF-alpha had a signifi- 
cant correlation among themselves. Likewise, MMP-9 
correlated significantly to VEGF, (r = 0.51) and inversely 
to inflammatory biomarkers. Primary chemoattractants, 
MCP-1 and PARC/CCL18 exhibited very little relation- 
ship to the other biomarkers, but their levels correlated 
weakly (r = 0.24). 



Biomarkers and clinical outcomes 

The differences in clinical outcomes between the 
patients in the high and low quartiles of biomarker levels 
are expressed using the OSCAR plot are shown in 
Figure 2. Mostly, minor differences were observed in 
demographics between subjects with high or low quar- 
tiles of biomarker levels. However, the patients in the 
upper quartile of inflammatory markers (IB) showed 
more severe clinical impairment (red color in the 
OSCAR plot). The difference in the mean value of 
physiologic and QOL measurements between the upper 
and lower quartiles shows a lower lung function (FEVi: 
-138 ml, DLCO: -14%) and walking distance (- 59 m), 
worse health status (10.7 pts), BODE score (1.75 
points), and survival (-21%). In contrast, the subjects in 
the upper quartile of markers related to injury and re- 
pair (IRB) seemed to have a less severe disease status 
(shown in green color): better FEVi: 345 ml, DLCO 
19%, 6MWT: 83.5 m, SGRQ: - 18.4 points, BODE -1.5 
points and improved survival (28.5%) These differences 
are within the recognized MCID for those outcomes 
(data analysis section). The association between PARC/ 
CCL18 and MCP-1 and disease status was less evident. 

Using C statistic the biomarkers by themselves pre- 
dicted mortality with a value of 0.78 (MMP-9 and IL-6 
having a significant contribution). Inclusion of the 
BODE index improved the model predictive value to 
0.85 (MMP-9 was the only marker that provided signifi- 
cant contribution to this model). 

Conclusion 

This study of a panel of serum biomarkers related to in- 
flammation, injury and repair, and chemoattractants, 
measured in a cohort of patients with COPD has two 
important findings. First, the biomarkers level relate to 
the degree of airflow obstruction, functional capacity 
and health status. Second, a biomarker pattern (inflam- 
mation and destruction and repair) distinctly correlates 
with meaningful clinical outcomes including mortality. 
The interaction between the pattern and the outcomes 
can be expressed using a novel heatmap method 
(OSCAR plot). Taken together, these results suggest that 
a pattern based on values of a small selected panel of 
serum biomarkers can provide insight into the clinical 
expression of COPD. 

The validity and usefulness of biomarkers depend on 
several characteristics including: association to the 
pathophysiological processes, association to important 
clinical outcomes and sensitivity to detect clinically im- 
portant differences. [4,6] This study fulfills all of these 
characteristics. 

The serum level of the markers did relate to the 
severity of COPD measured by the degree of airflow 
limitation, functional capacity and health status. This 
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Table 1 Baseline characteristics grouped by GOLD stage 


Parameter 


All 


GOLD 1 


GOLD II 


GOLD III 


GOLD IV 


P- 
value* 


N 


253 


18 


73 


105 


57 




Age (Y) 


55 (9) 


59(13) 


66 (9) 


67(8) 


63 (9) 


0.005 


Females (%) 


36 


56 


38 


35 


30 


0.25 


Current Smoker (%) 


69 


53 


36 


25 


23 


0.059 


Smoking (p/y), median [IQR] 


60 [41-80] 


41 [35-60] 


60 [45-85] 


60 [42-80] 


66 [48-53] 


0.111 


BMI (kg/m^) 


27 (5) 


26 (6) 


29 (5) 


27 (5) 


25 (5) 


0.0002 


Fat Free IVIass Index 


1 9.6 (3.4) 


NA 


194 (1.9) 


1 9.7 (3.6) 


1 9.5 (3.9) 




Charlson Index 


0.5 (0.8) 


04 (0.9) 


0.6 (0.9) 


04 (0.8) 


0.3 (0.6) 


0.20 


FEV, (L) 


1.19 (0.54) 


2.25 (0.53) 


1.54 (041) 


1 .03 (0.26) 


0.65 (0.17) 


<.0001 


FEV, (%) 


46 (20) 


91 (10) 


61 (9) 


39 (6) 


23 (4) 


<.0001 


FEV, /FVC (%) 


46 (1 2) 


63 (6) 


54 (9) 


44 (10) 


35 (7) 


<.0001 


TLC (%) 


121 (23) 


120 (13) 


1 1 3 (20) 


1 20 (24) 


132 (25) 


<.0001 


FRC (%) 


156 (44) 


1 26 (23) 


132 (32) 


158 (37) 


1 94 (46) 


<.0001 


RV (%) 


212 (64) 


138 (14) 


173 (55) 


1 97 (56) 


260 (54) 


<.0001 


IC/TLC 


0.30 (0.1) 


0.44 (0.09) 


0.36 (0.07) 


0.28 (0.06) 


0.19 (0.06) 


<.0001 


SGRQ (total), median [IQR] 


44 [22-59] 


10 [3-19] 


29 [16-42] 


49 [34-59] 


60 [45-71] 


<.0001 


mMRC 2 and higher (%) 


52 


11 


22 


65 


83 


<.0001 


P02 (blood) 


70.6 (1 3.4) 


83.1 (10.9) 


71.5 (11) 


70.3 (13.6) 


65.6 (14.2) 


<.0001 


6 MWT (m) 


422 (118) 


530 (81) 


475 (99) 


425 (97) 


314 (105) 


<.0001 


BODE, median [IQR] 


3 [1-5] 


0 [0-1] 


1 [0-3] 


3 [2-4] 


6 [4-7] 


<.0001 


Exacerbation (%), preceding yr. 


60 


25 


59 


64 


64 


0.025 


Death during follow-up (%) 


35 


0 


17 


41 


60 


<.0001 


Continuous systemic steroid use (%) 


8 


0 


1 


10 


16 


0.012 


iL-D ipg/mLy 


^ A n Q 1/1/11 
0.4 [z.o- 1 4.4J 


/I n rn /I i 7 f^i 
4.U [U.4- 1 / .dJ 


4.9 [0.4-15.3] 


5.8 [2.9-10.6] 


y.3 [D.y-z 1 .uj 


U.U 1 


II -R fnn/iYil ^ 

IL O Vp^/IIILy 




y.u [o. 1 1 o.yj 


9.1 [6.2-13.9] 


9.8 [6.6-14.7] 


1 1 .Z [/ . J ZU.H-J 


061 


IL-16 (pg/mL) 


331 [239-547] 


314 [272-359] 


312 [209-430] 


347 [233-594] 


385 [294-635] 


0.052 


TNF-alpha (pg/mL) 


18.1 [8.5-67.8] 


13.1 [7.4-19.3] 


1 1 .7 [7.6-32.4] 


19.3 [9.5-92.2] 


26.5 [10.6-1048] 


0.005 


MMP-9 (pg/mL) 


5767 [3057-20700] 1 


8088 [9137-39533] 


0311 [3531-24100] 4675 [2975-18200] 


4214 [2975-13533] 


0.003 


VEGF (pg/mL) 


71.1 [31.7-165.4] 


100.5 [79.8-165.1] 


1 1 7.2 [34.8-238.2] 


52.0 [30.6-154.6] 


42.1 [27.4-83.5] 


0.001 


PARC (pg/mL) 


52372 [34358- 


541 77 [35.640- 


51358 [41459- 


50036 [29334- 


54232 [34358- 


0.7204 




66350] 


62458] 


64272] 


67443] 


77318] 




MCP-1 (pg/mL) 


527 [410-711] 


567 [435-642] 


484 [380-586] 


547 [440-775] 


609 [377-762] 


0.0169 



Footnotes: 

Values are means (SD) unless otherwise stated. Biomarkers are expressed median [IQR]. 
* p-value comparison across GOLD stage categories. 
NA = measurements not available. 



relationship varied depending on the function of the 
biomarker. Whereas the serum levels of the inflamma- 
tory markers IL-6, IL-8, IL-16, TNF alpha were higher, 
the markers of injury and repair MMP-9 and VEGF 
were lower in patients with more severe disease 
(Table 1). Some of our findings are consistent with 
those reported in previous cross sectional studies. Bon 
and collaborators described a statistically significant as- 
sociation between the degree of emphysema by chest 
CT and serum level of IL-6 and TNF alpha in 234 



COPD subjects. [32] TNF alpha is also known to in- 
duce synthesis of IL-16 [33] a marker mostly asso- 
ciated with asthma but generated by bronchial 
epithelium and dendritic cells, CD8-I- and CD4-I- T 
cells, key elements in COPD pathobiology [28]. Serum 
levels of IL-8, IL-6 and TNF alpha decreased in 
patients with COPD after lung volume reduction sur- 
gery [34] and their decrease correlated with the reduc- 
tion in hyperinflation and improvement in body mass 
composition post surgery. 
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Distribution of IL-6, IL-8 and TNF-alpha (median) by GOLD stage In COPD 
patients during the visit 1 (values in pgAnL) 
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Figure 1 A-C. Distribution of biomarkers values according to GOLD stage. 



□ MMP-9 
■ PARC 



Likewise, the inverse relation between the level of VEGF 
and worse COPD severity follows previous reports. VEGF 
receptor blockage has been shown to generate emphy- 
sema suggesting a role of VEGF in preventing destruction 
or helping facilitate lung repair [35]. Likewise, Valipour 



et al. [36] found a significant correlation (r = 0.47 
p < 0.001) between the degree of airway obstruction and 
VEGF serum level in 30 stable COPD patients. 

On the other hand, the association of MMP-9 and 
lung function differs from previous studies with limited 



Table 2 Associations among the biomarkers at baseline 



Biomarkers 


IL-6 


IL-8 


IL-16 


MCP-1 


MM9-9 


PARC 


TNF-Alpha 


VEGF 


IL-6 


1 


0.53 


0.06 


0.11 


-0.06 


0.09 


0.34 


0.09 


IL-8 


0.53 


1 


0.16 


0.14 


0.04 


0.10 


0.35 


0.08 


IL-16 


0.06 


0.16 


1 


0.05 


-0.18 


-0.08 


0.52 


-0.35 


MCP-1 


0.11 


0.14 


0.05 


1 


0.05 


0.24 


0.12 


0.14 


IVlMP-9 


-0.06 


0.04 


-0.18 


0.05 


1 


0.11 


-0.22 


0.51 


PARC 


0.09 


0.10 


-0.08 


0.24 


0.11 


1 


-0.07 


0.19 


TNF-Alpha 


0.34 


0.35 


0.52 


0.12 


-0.22 


-0.07 


1 


-0.22 


VEGF 


0.09 


0.08 


-0.35 


0.14 


0.51 


0.19 


-0.22 


1 



Significant associations (p< 0.001) are denoted with bold characters. 
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Inflammation 


Repair & Injury 


Chemoattractants 






IL6 




L8 


IL-16 


TNF alpha 


MMP-9 


VEGF 


PARC 


MCP-1 


Variables 


Upper 
QrtI 
(n=63) 


Lower 
QrtI 
(n=63) 


Upper 
QrtI 

(n=63) 


Lower 
QrtI 
(n^63) 


Upper 
QrtI 


Lower 
QrtI 


Upper 
QrtI 


Lower 
QrtI 
(n=62) 


Upper 
QrtI 

(n=62) 


Lower 
QrtI 
(n=62) 


Upper 
QrtI 

(n=62) 


Lower 
QrtI 
(n-62) 


Upper 
QrtI 
(n-61] 


Lower 
□rtl 

(n=62) 


Upper 
QrtI 
(1=82) 


Lower 
QrtI 
(n=61) 


Age in yrs (mean, SD) 


65 
(9) 


64 
(10) 


64 
(9) 


65 
(9) 


66 
(9) 


67 
(8) 


67 
(8) 


63 
(11) 


65 
(9) 


66 
(9) 


65 
(10) 


66 
(9) 


(10) 


(9) 


65 
(8) 


66 
(9) 


Pack- Years (mean, 
SD) 


58 
(27) 


58 
(24) 


62 
(32) 


58 
(27) 


65 
(41) 


61 
(24) 


59 
(35) 


60 
(24) 


63 
(30) 


56 
(35) 


65 
(28) 


64 
(37) 


68 
(34) 


60 
(37) 


63 
(33) 


62 
(30) 


Current Smoker (%) 


35 


33 


32 


30 


20 


27 


23 


37 


39 


27 


48 


19 


28 


32 


21 


20 


Males (%) 


70 


65 


68 65 


49 


81 


58 


57 


65 


63 


74 


57 ^ 




57 


57 


60 


BMI (mean, SD) 


26 

(5) 


27 

(6) 


26 

(6) 


28 
(5) 


27 

{5} 


28 
(5) 


27 

(5) 


27 
(5) 


28 

(e) 


26 
(5) 


28 
(5) 


26 
(6) 


26 
(5) 


27 

(5) 


27 

(5) 


26 
(5) 


FEV1 (L) (mean, SD) 


1.21 
(0.54) 


1.36 
(0.57) 


1.22 
(0.57) 


1.31 
(0.59) 


1.04 
(0.39) 


1.21 
(0.44) 


1.08 
(0.44) 


1.22 
(0.57) 


1.37 
(0.66) 


1.07 
(0.43) 


1.39 
(0.59) 


1.0 
(0.46) 


1.09 
(0.47) 


1.17 
(0.58) 


1.05 
(0.41) 


1.22 
(0.57) 


DLCO % (mean, SD) 


68.8 
(30.0) 


73.9 
(31.2) 


64.6 
(28.5) 


67.9 
(29.9) 


52.5 
(19.9) 


72.1 
(32.4) 


51.0 
(21.4) 


69.6 
(29.4) 


72.1 

(26.9) 


57.4 
(25.8) 


76.9 


52.3 


58.7 
(27.3) 


53.8 
(20.5) 


55.3 
(21 .1) 


60.7 
(27) 


IC-TLC (mean, SD) 


0.30 
(0.09) 


0.33 
(0.09) 


0.29 
(0.1) 


0.32 
(0.1) 


0.27 
(0.09) 


0.32 
(0.08) 


0.27 
(0.09) 


0.32 
(0.1) 


0.33 
(0.09) 


0.28 
(0.1) 


0.34 
(0.08) 


0.26 
(0.1) 


0.28 
(0.1) 


0.29 
(0.1) 


0.28 
(0.8) 


0.29 
(0.1) 


6MWD (mean, SD) 


420 
(123) 


476 
(86) 


418 

(130) 


447 
(110) 


389 
(125) 


463 
(94) 


377 
(120) 


453 
(111) 


469 
(96) 


402 
(120) 


480 
(89) 


380 
(125) 


381 
(124) 


424 
(110) 


382 
(114) 


435 
(114) 


SGRQ Total Score 
(mean, SD) 


40 
(23) 


34 
(21) 


38 
(25) 


36 
(22) 


49 
(21) 


32 
(19) 


47 
(21) 


38 
(20) 


30 
(19) 


49 

(18) 


31 
(20) 


48 
(21) 


46 
(21) 


48 
(21) 


53 
(22) 


42 
(22) 


BODE [median, IQR] 


3 [1-51 


1 [0-4J 


3 [1-5] 


2 [1-4] 


4 [2-6] 


2 [1-4] 


4 [3-5] 


2 [1-3] 


2 [1-4] 


3 [2-6] 


2 [13] 


4 [2-5] 


3 [2-5] 


4 [2-5] 


4 [2-6] 


2.5 [1-5] 


Mortality (%) 






39 


31 


42 


32 






20 49 


23 


51 




42 


33 


Charlson CI > 0 (%) 


30 


33 


34 


26 


29 


33 


25 


26 


34 


37 


31 


32 


34 


30 


29 


36 



The difference between the quartlles was considered as significant based on known MClDs for clinical measures or a comparison using chl-square, t-test or wllcoxon tests (p<0.05) for 
demographics and medical history. 

Figure 2 Differences in clinical outcomes between the patients in the high and low quartlles of biomarker levels (Observed Clinical 
Association Results or OSCAR plot). 



phenotypic information. We previously reported a direct 
association between the level of MMP-9 and COPD 
exacerbations [7]. Higashimoto et al. [37] and Olafsdottir 
and co-workers [38] described a negative correlation be- 
tween the serum level of MMP-9 and FEVi (r = -0.28 
and r = -0.11 p < 0.01). However Bolton et al. [39], found 
no correlation between the level of MMP-9 with FEVi in 
70 patients with COPD but higher values in those 
patients with concomitant osteoporosis. Furthermore, a 
low ratio of MMP-9/TIMP-1 (tissue inhibitor of metallo- 
proteinase 1) has been associated with airway obstruc- 
tion, suggesting that an excess of TIMP-1 compared to 
MMP-9 could be responsible for airway remodeling [40]. 
Unfortunately, we did not measure the level of TIMP-1, 
a regulator of MMP-9 activity, but it is possible that a 
dysregulation in this balance could explain our results. 
Interestingly, the role of matrix metalloproteinase has 
been recently reviewed, suggesting a role of MMP-9 not 
only in alveolar injury but also on the alveolar epithelial 
repair process in acute lung injury [41]. Although some- 
what unexpected, the consistency of our observations 
and the findings by others [39-41] suggest that we may 



have to re-evaluate the current thinking about the func- 
tion of MMP-9 in COPD and suggests a need to study 
the ratio MMP-9/TIMP-1 in COPD as well as in the 
concomitant presence of osteoporosis. 

The second and perhaps most important finding in 
this study is the relationship between serum biomarker 
levels and the clinical manifestations of the disease. We 
developed a novel expression of this relationship using 
the OSCAR plot, where the differences in clinical im- 
portant outcomes are compared between patients with 
the highest biomarker levels and those with the lowest 
quartile. The differences are not only statistically signifi- 
cant but have recognized clinical relevance. The OSCAR 
plot provides a comprehensive visualization that relates 
biomarkers levels patterns to outcomes thus presenting 
a "fingerprint" of the biomarker/clinical relationship. 
The results suggest that the patients with the highest 
levels of inflammatory markers have more severe disease 
and higher mortality (concepts in line with current 
thinking). Patients with the highest level of VEGF and 
MMP-9 had better physiological and functional out- 
comes. We believe this type of analysis is useful in 
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providing patterns profiles that may help express the 
complexity of COPD. Even though we determined by 
multiple regression analysis that decreased levels of 
MMP-9 were consistently associated with the increased 
odd of mortality, we believe a pattern profile using 
OSCAR plots may be more informative than associa- 
tions between single biomarker level and individual out- 
come in COPD. 

This study had some limitations. First, we acknow- 
ledge that the selection of biomarkers is incomplete, but 
we believe is reflective of several mechanistic pathways 
including inflammation, chemoattraction and injury and 
repair. The selection of the markers is based on previous 
studies that included high throughput proteomic analysis 
of 147 markers in stable and unstable conditions and 
smokers and non-smokers controls [7,42]. The ECLIPSE 
study also reported differences between controls and 
COPD patients using 4 of the biomarkers here included 
(IL-6, IL-8, MMP-9 and CCL-18) [12] [43]. Second, al- 
though it could be argued that there was no validation 
cohort, this study follows our previous work in serum 
biomarkers but in a multicenter fashion, with a much 
larger set of patients and adding mortality as an out- 
come. Third, the relative low comorbidity rate of this co- 
hort of patients may not represent other COPD cohorts. 
However, they were recruited from 2 COPD centers in 2 
different countries, with similar COPD severity by 
GOLD and BODE index [43] and the results may there- 
fore reflects findings more likely associated to COPD 
than to other comorbidities. Fourth, there was no radio- 
logic evaluation (chest computed tomography) in this 
study to correlate changes in emphysema scores and air- 
way thickness with biomarkers level. However, we per- 
formed diffusion capacity measurements (DLCO), a 
marker of emphysema severity. The level of DLCO cor- 
related with several biomarkers as shown in Figure 2. 
Fifth, we acknowledge this to be an exploratory study 
and one that needs to be replicated in a larger cohort. 
The statistical analysis of the biomarker results by quar- 
tiles allowed us to determine clinically important differ- 
ences not otherwise capture by studying the mean value 
for the entire group. Since COPD is a heterogenous dis- 
ease, elevated biomarker level is not universal and there- 
fore, focusing on those individuals with the highest 
inflammatory level and/or the lowest repair biomarker 
level may be informative and of clinical importance. Fi- 
nally, the relatively small number of patients prevents a 
complete exploration of the value of individual serum 
markers to predict outcome above and beyond well 
known clinical predictors such as BODE. The relatively 
small cohort size may have also affected the predictive 
value of PARC/CCL18 that was found to be associated 
with mortality in two recently large cohorts (Lung 
Health Study and ECLIPSE with 4,825 and 1809 subjects 



respectively). However, the association between IL-6 and 
MMP-9 with three-year mortality in this study of a 
smaller size indicates a possible role for these biomar- 
kers in larger cohorts [44]. 

In conclusion, we have shown in a cohort of patients 
with COPD that the patterns of a panel of serum bio- 
markers levels relate to disease severity and outcome. 
The patterns support the role of inflammation as one 
driver of disease severity and a possible protective role 
for elevated markers of tissue repair. It may very well be 
that the balance between both processes determines the 
final outcome of individual patients. The novel way used 
to evaluate the relationship between serum biomarker 
levels and clinical outcomes (OSCAR plots) may help 
explore avenues of research related to the different 
dimensions of the disease. 
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